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Abstract 

Coincidence probabilities, which yield Renyi entropies, are inves¬ 
tigated in a generalized Gaussian model, which includes interparticle 
correlations. 
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1 Introduction 

Let us consider an M-particle state and assume that its density matrix is 

p{ K, q) = e -«(K)-|qL 2 q+iqX(K) j (l) 

where K = |(p + p’), q = p — p’ and X(K) are 3M-dimensional vectors, 
and L 2 is a 3 M x 3 M dimensional matrix. Matrix L 2 is in general a function 


1 


of K. The physical interpretation of formula Q is naturally obtained when 
the density matrix is converted into the corresponding Wigner function: 

/9T 3M — — 

W( K, X) = V g—f(K) —i(X—X(K))L~ 2 (X—X(K)) (2) 

DetL 

Note that DetL is an effective volume of the system. 

In this paper we derive, assuming (JTJ) , an asymptotic formula, valid when 
the eigenvalues of matrix L are large, for the coincidence probabilities 


<7(0 = Trp 1 


d l p fJp(K i ,q i 


3 = 1 


which are needed to obtain the Renyi entropies 


(3) 


Hi = jD iog c(l). (4) 

Here K, = |(p ? + p J+1 ), q ; = p^- - p j+1 and p l+1 = p,. 

As seen from formula (HD, when the eigenvalues of matrix L are large, 
the integral is dominated by the region q ~ 0. Therefore, we propose to 
expand the exponent in the integrand in powers of the components of q and 
to keep only the terms up to second order. Thus, the problem reduces to the 
evaluation of a Gaussian integral. 

Our final result is an elegant formula generalizing the formulae derived 
for the one-dimensional Gaussian model by a number of authors 0-0- For 
the convenience of the reader, a short derivation of the formula for the one¬ 
dimensional case is given in Section 3. In Section 2 the integral © is worked 
out. In Section 4 this integral is significantly simplified using the result 
derived in Section 3. Section 5 contains our conclusions. Two Appendices 
contain discussions of the matrices R 2 and S introduced in the text. 


2 Evaluation of C(l) 

In order to calculate the coincidence probabilities it is convenient to introduce 
the notation 


Kj = K + k,, 


K 


1 

7 


(Pi + --- + Pi). 


(5) 
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Note that 


i i 

5^ k j = °; 5^' = 0 - ( 6 ) 

3 =1 3 =1 

Aii immediate implication is that the terms with A" do not contribute. In¬ 
deed, let us consider the expansions of each term gjX(Kj) in powers of the 
components of q. The first terms are g.,X(K) and their total contribution to 
the exponent is 


i 

iX(K)£q, = 0. (7) 

3 = 1 

The total contribution from the second terms of the expansion must be zero 
because of the hermiticity of the density matrix and all the further terms are 
negligible being cubic or higher order in the components of q. 

It is convenient to replace in © the variables p l5 ..., by the variables 
K. q l5 ..., q^. The Jacobian of this change of variables equals one. Further, 
one introduces into the integrand the factor 



where t is another 3 M- dimensional vector, and compensates it by the inte¬ 
gration over q-. Thus, © takes the form 

C(l)=[dK [ T ^— . ( 9 ) 

J J (27 xyM J 

The sum of the v terms can be rewritten as follows 

i 1 i 

X>( k j) = MK) + ^ J>;V)(kWMK), (10) 

3 = 1 i =1 

where the differential operators V Q act on the components of K. Expressing 
the vectors k ? as linear combinations of the vectors q - and performing the 
summation over j one Ends that 

i 1 i -i 

V v(Kj) = lv( K) + - qjSjkR 2 qf., (11) 

3=1 j,k 
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where R 2 is a 3 M x 3 M matrix with elements 


li 2 a0 = V a V^(K) (12) 

and S is an l x l matrix. Matrix R 2 for a specihc model is discussed in 
Appendix 1. The matrix elements Sj k are not needed for our calculation. 
They are, however, useful for cross-checks. Therefore, they are calculated 
and discussed in Appendix 2. Actually, as seen there, there are potentially 
useful ambiguities in the definitions of the matrices S. 

The next step is to replace the integration variables q- by their linear 

combinations Q •, which are the eigenvectors of matrix S: q^ = Y^k=i UjkQk- 
This transformation can be chosen orthogonal and then its Jakobian equals 
one. Using the notation 


U T SU = A; A jk = Aj8 jk 

one finds 


(13) 


——Q(L“ + AjR“)Q + it 'y ^ UfcjQ 

k =1 

(14) 

Further, for each j we introduce an orthogonal transformation Q = V^z 
which diagonalizes the quadratic form in the exponent 

V^ t (L 2 + A,R 2 )UW = A^; A$ = A ^8 a0 . (15) 

This reduces each integral over Q into a product of 3 M single Gaussian 
integrals. Performing them one gets 


C{1)= dK e~ lv(K) 


dt 


(2tt 


13 M 


l 

n 

3 =1 


' d,Q exp 


C(f) = V2? M(M) JdK 
where 


e -lv( K) 


n ; =1 Det (L 2 + A ; R 2 


cj = 7 J2 u v u > 


kj 


ik 



2 ] A0) 

3 =1 

(16) 

(17) 
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and the identities 


3 M 

n>“> = I«(L 2 + A,R 2 ) (18) 

a=l 

have been used. 

Performing in m the integration over t and extracting DetL 2 = ( Deth ) 2 
one gets 


C(l) = 


where 


f 

[DetL J 


3M(l—l) 


l 


3M 
2 


dK 


0 -lv( K) 


DetA nhi Det( 1 + \ ; l. ‘R 2 L 


(19) 


A = ;^QL-va)^l 7 vW T L- 1 . 

3 = 1 

Using (usd this can be also rewritten as 

A = /V - ^"2-T 

3 = 1 J 

which yields 


( 20 ) 


( 21 ) 


C(i) = 



3M(!-1) 


l 


3M 
2 


7/K 


e —MK) 




EU^nf’d+A.x-'^L- 


( 22 ) 

The symbol nL=i denotes the product nL=i with the j-th factor omitted. 
This is a usable expression for the coincidence probabilities C(Z) 1 . The for¬ 
mula, however, can be significantly simplified. In order to derive the simplifi¬ 
cation we will need a result concerning the one-dimensional Gaussian model. 
This is derived in the following section. 


Tt requires, however, the determination of the eigenvalues A& and of the coefficients 
Cj. This is a nontrivial task for large l. 
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3 The one-dimensional Gaussian model 


The model considered in the present paper is a generalization of the more 
common Gaussian model, where the M-particle density matrix is a product 
of 3 M one dimensional Gaussian density matrices of the form: 


P(K,q) 


/; In-K- ii-./ 


(23) 


where K , q. R , L are just real numbers. Let us calculate the coincidence prob¬ 
abilities Ci{l) for the one-dimensional model (12311 . To this end we observe 
that p(K, q ) is diagonal in the representation of the wave functions of the 
harmonic oscillator 0-0 and can be written in the form 


where 


P(K, = V ; n(p)AnV’nb')> (24) 

n =0 


V>n(p) 


Therefore, 


= (1 - z)z n ; 


a 


-^exp (~l aV ] H n (apy, 


= VRL ; 


z = 


1 - 2A 

2L 

1 ' 2L 


Trp' = E = (1 " z )‘ 


71=0 


l-V 


(25) 

(26) 

(27) 


(28) 


where a geometric progression has been summed. It is seen from (E3) that 
|z| < 1, as it should. 

To show that formulae (1231) and (1231) are equivalent, one may use the 
identity 2 


E 

71=0 


z n 1 

—H n (x)H n (y) = -^== 


Az 2 


exp 


-4z 


z(x 2 + y 2 ) — xy 
1-4 z 2 


(29) 


2 To prove it is enough to substitute on the left hand side twice the definition 
H n (u ) = fpjlt (u + * , perform the summation and a Gaussian integration. 
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Substituting the formula for z into (EBP one easily finds 


Ci(l) = 

This is the formula we need 


R\ l 


L > (i+ &f-(i-£)‘ 


( 30 ) 


4 Final expression for C(l) 

The one-dimensional Gaussian model can be also studied using the methods 
from Section 2. There is a number of simplifications, however. Since the 
dimension of the matices 3 M is replaced by one 


A (j) = L 2 + AjR 2 - (d= 1). (31) 

Since there is no need to diagonalize L 2 + A*i? 2 : V^ J) = V^ 1 = 1. Moreover, 


e -iv( K ) = (JL V e -| 


r 2 k . 




(d= 1). 


(32) 


Tims 


Ci(l) 


R l 

toy/UU(L 2 + ^R 2 ) 


' (IK e 


-±i< 2 r 2 



2 ^ Cj L 2 + KR 2 
j =i J 


(33) 

Note that in spite of all these simplifications the eigenvalues A j and the co¬ 
efficients Cj are the same as in the general case. Performing the integrations 
over t and K\ 


where 


0,(1) = 


R 


• i -i 


AMnLm+A^ 2 ) 




c, 


L 2 + KR 21 

j =i J 


(d= !)• 


(34) 


(35) 
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The formula for C\{T) can be rewritten as 


Ci (0 = 


i?y _i i i 

L ' V^-iCjnSid+A^f) 2 )’ 


( 36 ) 


Let us compare now this result with formula dSHD- The equivalence of the 
two formulae implies that 


i Kj) 


l 


N 


£<?, n<i+v 



or equivalently 




(37) 


N 


1 1(3) 

E^rid+A) 


j =i fc=i 



N(l) 

1 + ^ 2 _2n 
n=l 


(/-!)! 


(2n + !)!(/ - 1 - 2n)! 
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(38) 

where iV(/) = E [|(l — 1)]. Note that this identity is valid whatever is 
substituted for -p-. 

Let us now go back to the general case. The argument of the determi¬ 
nant in formula differs from the argument of the square root in formula 
(Id(11) only by the substitution of the matrix L _1 R 2 L _1 for the number jj. 
Therefore, the same substitution can be made in the identity m and one 
obtains 


C(l) = 


where 


\ 

DetL ) 


\ 3M(I-1) 


I'T dK 


0 —lv[K) 


Det 


1 + E2 a n (L-'R'L- 1 ) 


1d2t -l\ n 


T, (39) 


Crt. — 2 


—2 n 


(/-!)! 


(2n+ !)!(/ - 1 - 2n)! 


(40) 


This is our final formula and the main result of this paper. Let us note that 
when the eigenvalues t±,..., t^M of matrix 


V = L' 1 R 2 L 1 
8 


( 41 ) 
















are known, the determinant in the integrand of formula m can be evaluated 
and the we get 


f e -MK>v / nKV 

DetL J 2 ./ dh [(1 + i (1 - i VQ‘]' 

5 Discussion 

We have derived an explicit formula for the coincidence probabilities in a 
general, multidimensional Gaussian model. This model can be considered 
as an approximation, at large volume of the system, of a model with an 
arbitrary momentum distribution. Therefore, our result may be useful for 
a rather wide class of physical situations, particularly for analyses of the 
systems created in heavy ion collisions. The formula obtained in the present 
paper is elegant, transparent and easy to use. 

Some comments are in order. 

• Since we do not assume factorization of the multiparticle density ma¬ 
trix, our calculation takes into account possible correlations between 
particles. Correlations show up as the non-diagonal terms in the ma¬ 
trix R 2 . There is no restriction on their character and magnitude. 

• It should be emphasized that, although the model is considered as an 
asymptotic expansion for large L, the terms in eq. El of higher order 
in L 2 should, in general, be included in spite of the fact that terms of 
these orders have been already neglected in the exponent of (JTJ). The 
point is that in the presence of correlations some eigenvalues of matrix 
L _1 R 2 L _1 may be of the order of the number of particles M. Then the 
expansion in formula (1HHI) contains terms of the order (^)”, which for 
high multiplicities are not necessarily small even for a large system. 

In order to illustrate the last point consider the matrices 3 

Rij = oi(5 i:j + /?); Lij = L5ij i, j = 1,..., 3 M. (43) 

3 The model could be made more realistic by distinguishing between the interparticle 
correlations and the correlations between the directions x,y,z. This, however, would not 
change the qualitative results we need. 
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For these matrices the eigenvalues of matrix L 2 R 2 L 1 are 
ai(l + 3 M/3) 


Ai = 


a 


L 2 


X2 — ■■■ — A3 m — ~j ~2 ' 


( 44 ) 


Therefore, 


Det 


W) 

l + J2 a n (L“ 1 R 2 L- 1 )’ 

n =1 


iv(0 




n=l 


a(l + 3M/3) \' 

a n | ^ ) 


N{1) 


1+E 


which shows that terms of the kind ( 3A ^ /3 )" do occur in the expansion. 


n=l 

(45) 


a 

n ^2n 


6 Appendix 1: Matrix V — L l R 2 L 1 in a 
model 


Let us consider a model where 


v(K) 

v(K) 
A 


M 




n=l 


Y 2 


-\/ K T + m 2 + — + log 
7T T(m + T)V8nA e~^. 


AE 


(46) 


In this formula T and A are parameters, m is the mass of the particle; K 
are the three-dimensional momentum vectors composing the 3 M dimensional 
vector K, K T is the transverse component of K\ the energy E = \/m'? + K 1 , 
the rapidity Y = A log ■ This is a commonly used model, where the par¬ 
ticles are uncorrelated, but there are some correlations between the x. y, z 
components of the momentum vectors. Thus, matrix R" consists of M diag¬ 
onal 3x3 blocks which will be denoted K 2 . The transverse momenta have 
a Boltzmann distribution, while in the longitudinal direction there is a mild 
cut-off by a Gaussian in rapidity. The last term in the expression for v en¬ 
sures the correct normalization. Matrix L is assumed to be diagonal with 
diagonal elements L x = L y = Lt, L z in every diagonal block corresponding to 
the diagonal blocks of R 2 . Further we use the labels a, j3 for the components 


3M—1 
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x,y,z and the labels a, b for the transverse components x,y, the transverse 
mass rriT = \Jm 2 + K\. 

Let us define 


This gives 


where 


V aP = 


R 


a/3 


V ab 

v az 

v zz 


L/yL 


a-L-ifi 


0 K a K h 

uK a , 

T, 


$5„ 


l 2 t g = 


1 2 

Tmj, E 4 

= 

1 1 
Tm T + E 2 

LtL z uj = 

1 (K z 

AE 2 \m 2 T 

L 2 ^ = 

1 

~A+1 

E 2 

A 


K z / K z 2 Y Y \ 

" AEm 2 T \Em 2 T + ^ + E*) ’ 
YK Z 


AEm 2 t ' 


Y\ 2 K z 
+ E E 4 


As seen from 


E EA 
-ED matrix T can be written in the form 


(47) 


(48) 

(49) 

(50) 


(51) 

(52) 

(53) 

(54) 


VaP ~ < f )| 5a/3 + VaP’ (55) 

The second row of matrix V' is proportional to the first. Therefore, DetV' = 
0. Consequently, one of the eigenvalues of matrix V equals $ and the other 
two can be found by solving a quadratic equation. The result is 


t\ = <h, (56) 

# 2 = \ (®K 2 + <f> + T + ^(QK 2 + $ — T) 2 + 4ca 2 A^ , (57) 

h = i f QKl + $ + T - sj(QK 2 + $ — T) 2 + Au 2 I\ 2 j . (58) 
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Since the eigenvalues are known formula © can be used to calculate C(l). 
In particular for l = 3 the result is 


C( 3) 



6 M 


3“ 


3 M 
2 



e -3v{K) 


/-, <t> W (eA'2+$)^-^2iC2 

+ I 2 M 1 + 12 


144 


(59) 


7 Appendix 2: Evaluation of the matrix S 


Let us represent each vector kj, j — 1,..., l as a linear combination of the 
vectors q^, j — 1, ,1 — 1 

1 1 

k j = 2 (Pj + Pj+i) “ y (Pi + ■ • • + Pz) = ( 60 ) 

k =1 

Because of the second identity © it is not mandatory to include qi in the 
expansion. Equating the coefficients of the vectors p l5 .. . on both sides of 
the second equality we get the equation system 


1 

2 


\&jn T ^j+l,n) 


1 

7 


Cjn 


n = 1,... ,i, 


(61) 


where c j0 = Cjj = 0 and 5i + \ >n stands for 5\ %n . Summing each side of the 
equations over n from n = k + 1 to n = l one gets from the resulting 
equalities 


c jk — —y + 1 + ~5ji — - [B(j — k) + @(j — k + 1)], 

where Q(n) equals one for n > 0 and zero otherwise. 

The coefficients S mn are dehned by 

i 

kjakj/3 ^ ^ S m nQmaQn0- 
m,n 


(62) 


(63) 


M ' 
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Substituting on the left hand side the expansions of k j in terms of the q y one 
finds 


Snn n 


E 

3 = 1 


CjmCjn- 


(64) 


Note that this matrix is symmetric and that by construction S m i = Si m = 0. 
Using the explicit formulae to perform the summation 


TTl(l — Tl) 1 

Smn = -7-“(1 + Smn), for Ul < U < l. (65) 

For m > n one finds the matrix elements from the symmetry S mn = S nrn . 

For l = 2 one gets 5 = 0, which greatly simplifies the analysis DUS- Note 
also that because of the second identity © one can add on the left hand side 
of relation (E2D a term 


i i 

E c b El CfcC fo 

j =0 k =1 


( 66 ) 


where Ck are arbitrary constants. This modifies S without affecting the prod¬ 
ucts kj a kjp and sometimes can be used to simplify the calculations. 

Choosing for instance: C\ = C2 = — C3 = — j-b for l = 3 and C\ = c 3 = 
—C 4 = |c 2 = | for l — 4 one gets the simple formulae: 


(t-2 0 o\ 

5(3) = 0 i 0 , 

Vo 0 j-J 


( \ 0 4 0 \ 

0 I 0 -I 

0 \ 0 

V 0 -I 0 I J 


(67) 


from which the eigenvalues are immediately visible. 
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